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Abstract 

We present a gauge model for the bimaximal mixing pattern among the 
neutrinos that explains both the atmospheric and solar neutrino data via large 
angle vacuum oscillation among the three known neutrinos. The model does 
not include righthanded neutrinos but additional Higgs triplets which acquire 
naturally small vev's due to the type II seesaw mechanism. A combination 
of global Le — Lf^ — Lj- and 5*3 symmetries constrain the mass matrix for 
both charged leptons and neutrinos in such a way that the bimaximal pattern 
emerges naturally at the tree level and needed splittings among neutrinos 
at the one loop level. This model predicts observable branching ratios for 
r — > fififj,, which could be used to test it. 

PACS: 14.60.Pq; ll:30.Hv; 12.15.Ff; 



I. INTRODUCTION 

The observation of a deficit as well as of a zenith angle dependence in the flux of atmo- 
spheric muon neutrinos by the Super-Kamiokande |jl]J^ collaboration has provided strong 
evidence that there are oscillations among the known neutrino species. The five solar neu- 
trino experiments have added to this sense of excitement by their long standing result 
that there is also a deficit of the solar neutrinos, which can be given a simple explanation in 
terms of neutrino oscillations 0. It thus appears certain that neutrinos have mass and they 
mix among each other. Although the details are fuzzy on the exact nature of the oscillations 
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needed for the purpose, several very interesting scenarios exist. In particular, if the only 
laboratory indication of the neutrino oscillation by the Los Alamos collaboration (LSND) 
is not included in the picture, there is a mixing scheme known as the bimaximal mixing, 
where both the solar and atmospheric data are explained by large mixing among the three 
known neutrinos In this picture, solar neutrino puzzle could either be solved via the 
large angle MSW mechanism p| or via the vacuum oscillation mechanism [^] depending on 
the mass difference between the muon and the electron neutrinos. In this paper we will as- 
sume the vacuum oscillation between the v^. and as the solution, which requires that their 
mass difference square must be ~ 10~^° eV^. The observed electron energy distribution as 
well as some hints of bi-annual variation of the solar neutrino flux by the Super-Kamiokande 
collaboration may be pointing in this direction. 

If we accept this particular resolution of the neutrino puzzles, two major theoretical 
challenges emerge: one, how does one get naturally a theory that leads to the bimaximal 
mixing matrix and two, how does the same framework explain the tiny mass difference square 
(~ 10^^° eV^) needed for the purpose without fine tuning of parameters ? Our goal in this 
paper is to provide a simple model that generates both these features of the neutrino physics. 
Note that having a neutrino mass matrix of the right form to generate the bimaximal mixing 
pattern is by itself not sufficient since the desired mixing matrix is a combination of both 
the neutrino and the charged lepton mixing matrices i.e., V\lJy. Often it is assumed that 
f/^ = 1 by appropriately choosing couplings in a theory to have certain values. This is of 
course not technically natural since radiative corrections could induce arbitrary values for 
those parameters thereby upsetting the neutrino mixing pattern. So what is really needed 
is (i) a neutrino Majorana mass matrix of the right form to generate the bimaximal mixing 
and (ii) a diagonal charged lepton mass matrix . It is the goal of this paper to present a 
model that satisfies both these criteria naturally. In this respect our model is different from 
others discussed in the literature (see later for detailed comparison). 

One of the key ingredient in our work is the type II seesaw mechanism where the vev of a 
triplet Higgs becomes ultrasmall due to the presence of a high scale in the theory The 



presence of additional global symmetries in the model lead to a pattern of neutrino masses 
that leads to the bimaximal mixing among neutrinos while keeping flavor mixing among the 
charged leptons to be zero so that the bimaximal pattern dictated by the neutrino mass 
matrix that emerges is indeed natural. 

Using the definition of the mixing matrix as 

^1 \ 

(1) 

the bimaximal mixing corresponds to the mixing matrix f/j, given by ^ 
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As far as the neutrino masses go, they may be fully or partially degenerate or hierarchical 
as long as the mass differences fit the desired values. As mentioned, a convincing theoretical 
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explanation of this elegant mixing pattern seems so far to have been elusive, although there 
exist many interesting attempts |11-15|. The problem becomes even more challenging when 



we demand that the solar neutrino puzzle be solved by the vacuum oscillation mechanism. 

In this letter, we use the type II seesaw mechanism in conjunction with the global sym- 
metry 5*3 X f/ (l)e_^_T- to show that both these properties can be realized in a natural manner. 

This leads us to a neutrino mass pattern where m^^ <^ m^^ 
eV and a generalized bimaximal pattern given by: 
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II. THE MODEL 



We consider an extension of the standard model where the fermion content is left un- 
altered but with a Higgs sector extended as follows: three doublets 0o,0i,02, two triplets 
with Y = 2 denoted by Ai^2 and a charged isosinglet with Y = +2. The model has an S3 
symmetry (i.e. permutation group on three elements), under which the particles are assigned 
as shown in Table I. 



Fields 


S3 transformation 




2 




2 


Le, Gr 


1 


(01,02) 


2 


00 


1 


(Al,A2) 


2 


7]^ 


1' 



Table I: Transformation properties of the fields in the model under the 6*3 group. 

We also impose an Lg — — symmetry on the model. The Yukawa part of the 
Lagrangian in the leptonic sector consistent with these symmetries can be written as: 



Cy = hi {l^fiR + L^Tr 



{Lf^HR - L-rTR)(f)i + (L^Tr + L-rHR)(f)2 



+ heLeeR(f)o + fLe {L^Ai + Lr A2) + fL^LrTj + h.c. 



(4) 



We will show later by a detailed examination of the Higgs potential for the system that 
there is a domain of parameters for which we get the following vevs of the fields: 



-'1,21 



,)=vo; (0?)=t;i; (0^) = 0. 



(5) 



Clearly, the pattern of (j) vevs leads to a diagonal mass matrix for the charged leptons whereas 
the A vev's leads to a Majorana mass for the neutrinos of the form: 



/ mi 1712 
mi 
\ m2 



(6) 
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As a consequence, diagonalization of the neutrino mass matrix is solely responsible for the 



neutrino mixings and one obtains the pattern given in the generalized bimaximal form |ll^JT5 
(see Eq. (|])), with tan 6* = 1712/ mi = f J/ff. 

Let us now address the question of neutrino masses. Clearly, to understand the small 
neutrino masses, one must have a tiny value for the vev of the A fields. This is achieved by 
the type II seesaw mechanism [|l^]. This is a generic mechanism, which can be illustrated 
by the following simple model that has only one (p and one A. Consider the following Higgs 
potential for this system |jTO|: 

V{(j), A) = M^A+A - fi^(f)^(f) + X4(j)^<j)f + AaA+A 

+A0AA"fA0V + ^a<a<^AV0 + h.c. (7) 

Let us choose ~ 100 GeV and M ~ M/^^^^ ^ /i; in this case, the vev of {(jP) ~ /i whereas 
the vev of (A'^) ~ ^ ^ /i. This mechanism has been labeled type II seesaw and we see that 
if M ~ lO^'' GeV, then we get (A°) ~ 0.14 eV. In the presence of more A fields and extra 
symmetries that our model has this mechanism still operates and we have a small mass (in 
the 0.1 eV range) for the z/g and u^. The third eigenstate has zero mass. The z/^ and get 
mixed in the tree level and the mixing angle is near maximal unless we do fine tuning. As 
far as the and Ve are concerned, they are degenerate and have opposite CP; therefore at 
the tree level their mass difference vanishes. We will show below that their mass difference 
arises at the one loop level but due to the presence of the high mass scale that gave rise 
to the type II seesaw mechanism, the mass difference Amg_^ is naturally suppressed to the 
level of 10~^° eV^ without any unnatural fine tuning. The tree level mass matrix already 
explains the atmospheric neutrino puzzle due to the type II seesaw. 

Let us now turn to the explanation of the one loop contribution to the neutrino mass 
matrix. This is where the role of the 77 boson becomes important. Let us first note that the 
masses of the doublet bosons are of order of the electroweak scale ( 100 GeV) since their 
vevs must be of that order whereas that of the singlet rj and the triplet bosons are heavy 
(i.e., of order 10^^ GeV). In order to generate the mass difference between the z/g and z/^, we 
need nonzero entries for the fifi or ee element. Both of them will violate the Lf. — L^ — 
symmetry. This breaking is introduced by a soft term in the potential ?7*0i02 since 7] has 
Lf, — — Lt- = —2 and 0's are neutral under this global symmetry. This is a dimensional 
coupling and in accordance with our principle above that all fields which are not involved 
in the process of electroweak symmetry breaking are superheavy, we will choose this to be 
of order M. This leads to one loop graphs as in Figs. 1 and 2, which produce a neutrino 
mass matrix as follows: 



/ mi m2 
mi m^^ m^r I ; (8) 
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where m^^ ~ m^^ sin 9 cos 9 and 



_ fh2 mrM^^^v 

For M ~ M^0<^ and f'h2 ~ 10""^, we get m^^ ~ 10"^ eV, leading to Am^^ ~ 10~^° eV^, as 
is required to explain the solar neutrino puzzle via vacuum oscillation. Note that the tau 
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neutrino picks up very tiny mass at one loop (~ 10"^" eV). This completes the derivation 
of the main result of our paper. The important point to note is that we need to choose the 
Yukawa couplings /' and /12 individually only of order 10~^. 

Let us now compare our model with two existing ones [|12[|13| in the literature. The 
model of Ma [|T2| has a similar field content to ours in all respects except that there is only 
one triplet field as against two in our paper. Thus inspite of the S3 symmetry in that paper, 
the tree level mass matrix is very different and one does not have the bimaximal pattern 
at the tree level from the neutrino sector sector alone as in our case. On the other hand, 
the work of Ref. |T3[ uses only the — — Lj. symmetry, which allows arbitrary mixing 
angles in the charged lepton sector, which have to be set to zero at tree level. The addition 
of the S3 symmetry as we do in this paper helps us to keep the charged leptons diagonal. 
There are also other major differences between the work of and this work in the way 
the detailed dominant entries of the neutrino mass matrix arise- in our case at the tree level 
where as in |]T3 at the one loop model via the Zee-mechanism. 



III. RARE TAU DECAYS 

In this section we present a test of the model that involves flavor changing decays of 
the tau lepton. There are three sources of flavor changing effects in this model, i.e., via 
the exchanges of r/, A's and 02. At tree level, r/"*" and A^l^ exchanges lead to or r decay 
processes that include neutrinos as final products. However, as the masses of the exchange 
fields are very heavy (~ M) such processes are highly suppressed. Therefore, only the 
diagrams that involve the exchange of the doublet fields are important. 

The general exchange tree level diagrams for Flavor Changing Neutral Currents 
(FCNC) are depicted in Figures 3 and 4. The internal lines in those figures represent 
the contributions of the general mixings in the scalar sector that come from the trilinear 
and quartic scalar couplings. In both of the figures, since one scalar field is necessarily a 
02 field, the 5*3 symmetry restricts the involved vevs. In Fig. 3, for instance, it allows only 
(A2), regardless of whether the other field is 0o or 0i. Thus, we estimate the amplitude of 
this diagram to be 

hh,^-^^ ^ h . GeV-^; (10) 

where we have assumed as before that M/\^^ ~ M, ~ v, M^^^{A) ~ and taking 
^2 = TTiT-fv ^ 10~^. Here, h represent the Yukawa coupling of the vertex on the right. 

For the diagram in Fig. 4, the vevs could be either (0) or (A). Nevertheless, if we choose 
(0), again the 5*3 symmetry play an important role by constraining one of them to be (02), 
which is zero, then the only contribution comes from the (A) sector which is more suppressed 
already than the previous case by an extra 10~^^ times the quartic coupling constant, which 
being dimensionless may be chosen of the order of one. There is no compensating factor 
such as large M as in diagram in Fig. 3. 

Thus the only observable contribution to leptonic FCNC processes come from the di- 
agram in Fig. 3. Matching the external leptonic legs in Fig. 3 with the terms in the 
Lagrangian, we see that only observable processes are of type r — + /iyU/i and r fiee, which 
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have an amplitude of order hlO ^ GeV Summing up the two contributions and assuming 
(0o) — and all dimensional couplings to be same, we get h ~ ~ 10~^. This leads 

to the branching ratio for the decay mode B{t — > /i/x/x) — 10~^. Since in our estimate we 
have assumed several couplings to be of order one, the prediction is uncertain within an 
order of magnitude but we do not expect it to be much smaller. This may be compared with 
the present experimental bounds [|T^ for the decays r —>■ fififi which is 10^^. The branching 
ratio for the other allowed decay mode in our model B{t —>■ fiee) ~ 10^^^ is small due to 
the fact that it involves the electron Yukawa coupling which is ~ irie/mw- Therefore, the 
three /i rare mode is the only observable FCNC processes in r decays in this model and 
could therefore be used test. 

Let us stress that, an interesting feature of the present model is that the really rare 
process yU eee is automatically suppressed as it does not appear at the tree level. This is 
because, the only tree level coupling among electron and muon (or tau) involves A"*"^, which 
does not mix with 00, as would be required to get three electrons in the final states. 

At one loop order the most interesting process again appears in r physics, i.e. the rare 
decay t ^ fi'-f. The coresponding diagram is showed in Fig. 5. Now the decay width for 
this process is roughly estimated to be 



167r2 m\^rnl^^^ 



10-'^ GeV; (11) 



giving a branching ratio of about ~ 10 ^, which is again below the current experimental 
bound fT^ of 10"^ 



IV. ANALYSIS OF THE HIGGS POTENTIAL 



Lets us now show how (02) = arises naturally from the potential. Given the irreducible 
representations (irreps) of S^: 2^ = (xi,X2) and 2y = (^1,^/2), we may build the following 
singlets Ixy = xiyi + X2y2', ^'xy = ^i?/2 — X2yi', and the new doublet 2xy = {xiyi —X2y2, Xiy2 + 
X2yi)- Using this simple rules it is straightforward to find all possible 5*3 and gauge invariant 
terms that include the scalar fields of the model. Such potential may be decomposed as 
V = V{(f)) + V{A) + \^(0, A) + \^(0, A,?7). The last term involves all the expressions 
containing 77. They do not contribute to the minimization of the potential, thus, it is not 
necessary to show them explicitly. 

As we already discussed above, the type II seesaw formula arises from V{A) + V{(f),A) 
by assuming large trilinear couplings. In this case (A) becomes much much smaller than (0) 
and then we may neglect those terms in the analysis of the vevs. Defining $ = (0i,02), 
the relevant terms of the potential are then represented by 

V{(P) = I? + III 1^2 + Al (1$2)2 + A2 (1^2)2 + 

A3 1$2 l</,2 + A3 1$4 + A5(l'<j,2)2 + Ae 1$3 l^;,; (12) 

where 1^^ means the singlet built by using n x irreps. By examing V{(f)) we may see that 
all the terms except the last one obey an accidental U{l)a symmetry, which makes itself 
evident if we parametrize the vevs as 
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COS a 
sin a 



(13) 



Then, in terms of these parameters, the potential reduces to 

V{ct>) = V{v, (0o)) + A6(0o)^''cosa. (14) 

If the last term was absent, the minimum would exhibit a flat direction and an ensuing 
Goldstone boson. However, the last term breaks this extra symmetry explicitly and removes 
the flatness. Moreover, from this expression it is now straightforward to see that the potential 
gets its minimum for a = (tt) if Ae < (> 0), which means that (02) = 0, as we expected. 
It is worth pointing out that this special effect in the potential does not appear on the A 
sector, since it is just a consequence of the presence of the extra singlet 0o- As a matter 
of fact, V^(A) is totally U{1)q symmetric while V^(0, A) contains several terms that break 
explicitly such symmetry in a less trivial way, then avoiding a null value of (A2), and giving 
the pattern of neutrino masses. 



V. CONCLUSIONS 

We have presented a model of neutrino masses based on the permutation symmetry 6*3 in 
conjunction with the Le — L^ — Lr symmetry which leads to the bimaximal mixing pattern at 
the tree level and explain atmospheric oscillations through the hierarchy m^,^ -C mj,^ ~ m^^ . 
The naturalness of the bimaximal mixing follows from the fact that the tree level mass 
matrix of the charged leptons is diagonal while the neutrino Majorana mass has a specific 
form dictated by the type II seesaw mechanism and the above symmetries. The soft breaking 
of the Lg—Lf^ — Lr symmetry by the scalar potential through a coupling of the scalar doublets 
with an 5*3 odd charged scalar, rj, leads to the small entries in the neutrino mass matrix 
via radiative corrections. They are responsible for the small splitting in Am^^ needed to 
explain the solar neutrino deficit via vacuum oscillations. This model can be tested via the 
rare decay r fififi whose branching ratio is predicted to be not too far below the current 



experimental limit |17 
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FIGURES 
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FIG. 1. One loop correction that generates the diagonal mass term m^^. A similar diagram 
provides m^-T-. 
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FIG. 2. One loop correction that generates the mass term m 
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FIG. 3. Generical Feynman diagram responsible for FCC in the model at tree level. The 
external lines are leptonic fields. The internal lines could in principle be any one of the scalars 
involved in the trilinear couplings allowed by the 5*3 symmetry. 
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FIG. 4. Tree level FCC diagram involving the mixing produced by quartic scalar couplings. 
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FIG. 5. One loop diagram that produce the t ^ jjL'y decay in the model. 
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